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Several accidents in towing operations of barges or 
disabled ships in restricted and open waters have made 
necessary the investigation of the course keeping stability 
of towed vessels. In this work a non-linear mathematical 
model is used to simulate the slow surge, sway, and yaw 
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effect of geometric parameters of the system on the 
stability of equilibrium configurations is analyzed. 

It is shown that for certain choices of towing system 
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results in qualitatively different asymptotic response. The 
results of this study identify regions in the parameter 
space that lead to either safe operations or hazardous 


system response. 
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l. INTRODUCTION 


A. BACKGROUND 

A long history of towing accidents resulting in loss of 
life, damage to property, and pollution of the environment 
have prompted many studies into the dynamics of towing 
operations. Of primary concern was the motions of the towed 
vessel in the horizontal plane (yaw, sway, and surge). 
Excessive and unstable motions could lead to collisions and 
capsizing. The ability to predict the motion of a 
particular towing system would be of particular benefit to 
ship designers and towing operators, by identifying those 
Situations where the towing operations would be the safest, 
or those which must be avoided. 

Previous studies at the University of Michigan and 
elsewhere had developed mathematical models and numerical 
technigues to analyze towing dynamics, and had identified 
those parameters which are of primary importance to the 
stability of the towing system. The linear model usually 
used to describe ship motions [Ref. 1, Chapter 7] is 
inadequate for the towing problem. Non-linear models, as in 
(Ref. 2}, must be utilized to accurately describe the 


towing system. These studies had identified the position of 


the towline attachment point on the towed vessel and the 
towline tension as the most significant (and controllable) 


parameters of the towing system. 


B. PROBLEM CONDITIONS 

In this study, computer programs developed in [Ref. 3] 
were used to analyze the effect of different parameter 
combinations on the towed stability of three vessels. These 
programs use bifurcation to identify the unstable and 
stable regions of the parameter space. The principal 
parameters studied were (Fig.1): 


1. longitudinal position of the towline attachment point 
forward of the towed vessel’s center of gravity, xp; 


2. athwartships position of the towline attachment point 
port or starboard of the towed vessel’s centerline, 


YP; 
3. length of the towline, Lw. 

In the model used in this study, unlike [Ref. 2], the 
towline is modeled as an inextensible caternary, thus 
making towline tension a function of its length. The model 
conditions were: 

1. speed of towing vessel of 2 knots; 
2. towing vessel on steady course; 


3. calm seas, no wind; i.e., no external environmental 
forces. 


Characteristics of the towed vessel were inputted into the 
programs from a data file containing hydrodynamic 


coefficients, resistance data, towline characteristics, and 






MOWED YESSEL 


TOWING VESSEL 


: Figure 1. Problem Geometry 


skeg data, if applicable. The effect of asymmetrical forces 
acting on the towed vessel, such as the presence of a 
propellor or an environmental force, are introduced through 
a bias in the data file. All dimensions are 
nondimensionalized with respect to the towed ship’s length 
between perpendiculars (LBP). 

Three vessels were studied (Fig.2): 


1. a 191 foot barge with a skeg, with no propellor 
(1.e., no bias); 


2. a 1066 foot tanker with no skeg, but with a propellor 
(1.e., with a bias); 


3. the same barge as in 1), but without the skeg and 
with a propellor (i.e., with a bias). 


Unlike previous studies, this work includes the effect 
of athwartship position of the towline attachment point in 
the stability of the towing system. 

Chapter II provides background into the problem 
formulation and stability analysis used in this study. 
Chapter III presents the results of the analysis and 
discusses some practical aspects of these results. Chapter 
IV discusses the conclusions which can be made from the 
results of this study on the stability of the towing system 


and the use of the techniques used herein. 





¢ Barge with skeg, no bias 


* Barge with no skeg, with bias 
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* Tanker with bias 





Figure 2. Body Plans of Vessels Studied 


li. PROBLEM FORMULATION AND METHOD OF APPROACH 


Slow motions of a towed vessel in the horizontal plane 
are described by a system of six nonlinear, coupled, 
differential equations. [{Ref. 4 and 5] In its standard form 
this system is 


e 1 
= 
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In the above equations, x, denotes the sway velocity in 
surge (longitudinal motion) of the towed vessel, x5 the 
velocity in sway (lateral motion), x, the angular velocity 
in yaw (turning motion about the vertical axis), x, and x, 
the coordinates of the center of gravity of the towed vessel 
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with respect to an (x,y)-coordinate system moving with the 
towing vessel, and x¢ the towed vessel yaw angle. Further, U 
is the steady towing vessel velocity in the x-direction, xp 
and yp are the coordinates of the towline connection point 
on the towed vessel with respect to an (X,Y)-coordinate 
system with its origin at the towed vessel center of 
gravity, and R,, Ry are the towline restoring forces. The 
towing system configuration and notation conventions are 
shown in Figure 1. Expressions for F,, F5, F3 are derived by 
Taylor expansion in terms of the relative velocities xj, X5, 
x3 Of the towed vessel with respect to the water. In 
nonlinear analysis terms up to third order are used whereas 
terms beyond third order and second- and higher-order 
acceleration terms are uSually neglected. Subscripts u, v, r 
indicate derivative of force-moment component with respect 
to X,, X5, X3 respectively, and subscript c indicates 
propellor dependent terms., which represent a source of 
system asymmetry. These terms are zero in the absence of a 
propellor. Terms Xobc: Yabo: Naber Where a, b, c are dummy 
independent variables representing u, v, r, are usually 
called slow motion derivatives. In unsteady reference 
motion, slow-motion derivatives are considered as functions 
of the frequency of motion. In our study of slowly varying 
reference motions, we assume that slow motion derivatives 
are time independent. This is a good approximation for ships 
with usual hull shapes and moderate speeds. 
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R, and R,, denote restoring forces from the towline, and 


x y 
for a quasistatic towline response they are expressed as 
implicit functions of x,, x5, X¢- In this study, the model 
used for the towline is that of an inextensible heavy 
caternary with nonlinear force-displacement characteristics 
aS given in [Ref. 3]. 

In compact notation the above system of six ordinary 


differential equations is denoted as 
x=f (x) Gu) 


where x and f are six dimenSional vectors. To analyze the 
stability properties of (1), the first step is to identify 
the equilibrium configuration of the system. For this we 
have to solve a system of six nonlinear, coupled algebraic 


equations 
f(x) =0 (2) 


where x denotes an equilibrium configuration. It can be 
shown [Ref. 5] that system (2) has at most three solutions 
in x corresponding to three distinct equilibrium positions. 
In this study we concentrated our efforts on one of these 
equilibrium positions, namely the one which, in the absence 
of a bias in the system, corresponds to the towed vessel 
being located directly astern of the tow-tug. This is the 
most interesting in applications. Having computed x, its 
Stability properties can be established as follows: 
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Linearization of (1) around x leads to the linear system 

2=Az (3) 
where z represents the excursion from the equilibrium x, and 
A is a constant 6x6 matrix. If all eigenvalues of A have 
negative real parts, then x is stable; otherwise it is 
unstable. 

In this study, we performed parametric analysis of the 
central equilibrium in therm of towline length Lw, and the 
towing point coordinates with respect to the center of 
gravity of the towed vessel, xp and yp. These parameters can 
be easily changed before or during towing operations and can 
provide a means of passive control of the towing system. 
Parameter Lw directly affects the amount of tension 
developed by the towine. Parameter xp is directly related to 
the towline restoring force and moment. A small value for xp 
may not be able to provide adequate restoring moment and may 
not guarentee system stability. On the other hand a very 
large value for xp may result in over-compensation and 
therefore instability. Nonzero yp values result in a source 
of asymmetry introduced in the system. For a biased system 
(for example due to the presence of a propellor or 
environmental forces), it should be expected that an extra 
appropriate bias introduced via a nonzero yp helps 
counteract the effect of the former bias, and hence improve 


stability. 


EO 


The particular equilibrium position will lose its 
stability when an eigenvalue of the A matrix in (3) changes 
its sign from negative real part to positive real part. The 
case when a real eigenvalue crosses zero has ben analyzed in 
detail in [Ref. 5]. This corresponds to a static loss of 
stability with generation of additional equilibrium 
positions in the form of solution branching. In this study 
we focussed our attention on the case when a complex 
conjugate pair crosses the imaginary axis. This corresponds 
to a Hopf bifurcation: the particular equilibrium 
experiences a dynamic loss of stability and the system 
begins to oscillate. The resulting periodic solutions can be 
stable or unstable, but at any rate, such a situation is 


hazardous and should be avoided during towing operations. 
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lll. RESULTS AND DISCUSSION 


A. BARGE WITH SKEG 

The first vessel to be studied was an unpowered barge 
with a skeg aft. Since there is no propellor to introduce 
a bias, the barge has athwartship symmetry. 

1. Figure 3: Critical Real Part vs. xp 

Program TOWBIF1 calculates eigenvalues for specific 
Lw and yp, creates a file for each of six real and six 
imaginary parts, and creates a separate file containing the 
largest real part. The real part with the largest value is 
the critical indicator of the system’s stability: if it is 
greater than zero, the system will be unstable; if less 
than zero, the system will be stable. 

Figure 3 shows plots for the critical parts for 
yp=0.05 and three values for Lw. The region where the plot 
is greater than zero indicate that range of xp where the 
system is unstable. For example, for Lw=0.5 the critical 
real part is greater than zero for the range of xp=0.18 to 
xp=0.48, so the system is unstable within this range. 

Note that as Lw increases, the unstable range 
becomes smaller, until for Lw=3.0 there is no region 
greater than zero. Therefore, the system will be stable 
for all values of xp; i.e, the barge should exhibit no 
unstable motion. 
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2. Figure 4: yp vs. xp, Lw as Parameter 

Program TOWBIF2 does the same calculations as 
TOWBIF1 over a range of values of yp with a given Lw, 
instead of a single value of yp and Lw. In essence, Figure 
3 represents a cut of Figure 4 at a single value of yp and 
Lw. Unlike TOWBIF1, TOWBIF2 writes a point only where the 
critical real part changes sign. When plotted, these 
produce a curve delineating stable and unstable regions of 
the parameter space. This is the point of the process; we 
are more interested in finding what parameters produce 
stable or unstable system than the actual results of the 
equations of motion. 

Recalling Figure 3, The area inside the curves 
represent the unstable region. For example, for Lw=0.5, 
the system is unstable for all values of yp within the 
range xp=0.2 to xp=0.5. Increasing Lw first decreases the 
unstable range of xp for high yp, then decreases the 
unstable range of yp. For large Lw (>4.0), the unstable 
range virtually disappears. 

3. Figure 5: Lw vs. xp, yp as parameter 

Program TOWBIF3 performs the same calculations as 
TOWBIF2, but with Lw as the ordinate and yp as the 
parameter. Thus Figure 5 provides the same information as 


Figure 4 but with a different perspective. 
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As in Figure 4, the unstable region is inside the 
curves. It clearly shows how increasing Lw decreases 
the unstable range for a constant yp, as was evident in 
Figure 4. It also shows that, for constant Lw greater than 
about 0.7, increasing yp also decreases the unstable 
region. In the narrow range of Lw from 0.2 to 0.7, 
increaSing yp increases the extent of the unstable range of 
xp. This effect is apparent in Figure 4, but more 
dramatically presented in Figure 5. It would appear that 
using two views of the data would emphasize aspects of the 
curve that may be overlooked with one view. 

Since positive values of yp represent port side 
placement of the towline attachment point, and negative 
values starboard side placement, both positive and negative 
values for yp were studied. As expected from the port- 
starboard symmetry of the barge, curves for positive and 
negative values of yp were identical, and only positive 
values were presented here. 

From an operational point of view, one may conclude 
from these curves that for the unpowered barge, placing the 
towline on an attachment point to either side, as far 
forward as possible, will make the tow stable for the 
greatest range of towline length, but the towline should be 


kept no shorter than the length of the barge. 
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B. TANKER 
The second vessel studied was a tanker typical of those 
now in service. The effect of the tanker’s propellor makes 
the hull asymmetrical; this effect is represented by a bias 
included in the tanker data file. 
1. Figure 6: Critical Real part vs xp 

Figure 6 plots data generated from TOWBIF1 with 
yp=0.10 and two values for Lw. These plots show the stable 
region to be between the two zero values for the curve. 
Note that the stable region becomes smaller with increasing 
Lw. 

Note also that both curves are discontinuous in their 
Slopes. The critical real parts file 1s a composite of 
several results files, each of which is critical over a 
certain range. Each results file forms a smooth curve; 
thus the curves plotted on each TOWBIF1 figure may be 
combinations of the critical section of several results 
files. 

Finally, note that the stable region occurs over a 
Narrow range of xp, unlike the barge with skeg discussed 
earlier. 

2. Figure 7: Lw vs. xp, yp as Parameter 

Figure 7 was produced from data generated by TOWBIF3 

for positive values of yp. As was shown in Figure 6, the 


Stable region is inside the curves. The vertical line at 
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xp=0.404 iS a common crossing point for all curves. Each 
curve is formed by two cusps, with the upper cusp 
dominating with decreasing yp. Each cusp is the plot of 
different critical pair of eigenvalues; the "nose" in the 
curves is the point where they intersect. 

Note how the stable region gets smaller with 
decreasing yp for Lw less than 0.7, for example, with 
yp=0.0 and Lw greater than 1.0, there iS a very narrow 
range of xp where stability can be assured. 

3. Figure 8: Critical Real Part vs. xp 

TOWBIF1 was again used to form Figure 8, this time 
with one value for Lw (Lw=0.6) and three negative values 
for yp. The negative yp curves pass from stable to unstable 
regions, with the stable ranges for xp getting smaller as 
xp becomes more negative. 

As in Figure 6, the curves are composites of those 
results curves which are critical over a particular range 
SrexD . 

4. Figure 9: Lw vs. xp, Negative Values of yp as 
Parameters. 

Figure 9 data was generated from TOWBIF3, with yp=0.0 
curve included to provide continuity with Figure 7. 

The stable region gets smaller as yp decreases from 
0.0. At yp=-0.10 the "nose" between upper and lower cusps 


appears to be tipping up, with the region inside the "nose" 
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being stable. The curves plotted in Figure 8 were formed 
using a value of Lw which cut through this nose, thus 
forming the sinuous curves which pass in and out of the 
stable region. Note that for the most negative values of 
yp, the cusps have disappeared, and the stable range of xp 
is slightly increasing. 

5. Figure 10: Lw vs xp, Negative Values of yp as 
Parameters. 

Figure 10 is a "close-up" of Figure 9, focusing on 
what is happening around yp=-0.10. The lower cusp tips up 
and merges with the upper to form a single curve. Note how 
the stable range of xp virtually disappears for Lw greater 
than 0.5 for yp=-0.10 and -0.11. As was seen in Figure 9, 
stable range for xp for Lw greater than 0.5 reappears with 


yp less than -0.15. 


C. BARGE WITHOUT SKEG 

The third vessel was a self-propelled version of the 
barge studied in Section A (not under power during tow), 
but without the skeg. As with the tanker, the presence of 
the propellor, simulated by a bias in the data file, 
introduces port-starboard asymmetry. 

1. Figure 11: Critical Real Part vs xp for Lw=1.5 

Figure 11 shows curves for three values of yp 

(greater than zero, zero, and less than zero) and one value 


of Lw. The curves show the stable range of xp to be between 
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the zero crossing points of the critical real parts, as in 
the tanker case. Also similarly to the tanker, the stable 
region increases with increasing yp. These results are 


opposite to the propellor-less barge with the skeg. 


2. Figure 12: Lw vs xp 

Figure 12 dramatically shows how decreasing yp 
reduces the stable region. The vertical line at xp=0.14 was 
common to all vaiues of yp greater than and equal to Zero. 
For values of yp less than zero, the smooth shape of the 
curve is apparent. 

The tanker and self-propelled barge cases 
dramatically demonstrate the effect that a bias, like a 


propellor, can introduce to the stability of the system. 


D. PRACTICAL OBSERVATIONS 

Analysis of the graphs suggests some general principles 
which may be applied when conducting slow speed towing 
operations with the vessels discussed in this chapter. 
While these principles are of course not generally 
applicable to all vessels, they illustrate how the analysis 
techniques employed in this work can be applied to other 
vessels. 

For the unpowered, symmetric barge with a skeg, the 
operator should have the towline attachment point as far 


out to either side as possible and the towline as long as 
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practical. The attachment point can then be placed at any 
location forward of the center of gravity with stability 
assured. Conversely, if the attachment point must be on 
the centerline, placing it as far forward as possible 
(about half the barge’s length forward of the center of 
gravity) will assure stability for all towline lengths. 

For vessels with an asymmetrical bias (e.g., with a 
propellor), but without skegs, the attachment point needs 
to be as far to the biased side as possible (in the cases 
of the tanker and self-propelled barge, the +yp or port 
Side) and placed forward of the center of gravity the 
distance indicated on the graph for all towline lengths. 
Placement of the attachment point on the opposite side (in 
the cases studied, the starboard side) will virtually 


assure the system to be unstable. 


29 


IV. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

This study highlighted the effect of athwartship 
position of the towline attachment point. The common 
assumption among ship operators prior to this reseach held 
that placing the towline on the centerline on the foremost 
point of the towed vessel would create the most 
advantageous towing situation. Studies such as [Ref. 2] 
have shown that towing stability can be dependent on the 
longitudinal placement of the attachment point. This 
research has shown that for certain conditions, attaching 
the towline off the centerline can also improve towed 
stability. The optimum towing configuration requires a 
combination of all three parameters - longitudinal and 
athwartship placement of the towline attachment point, and 
towline length. 

The bifurcation technique used in this study can be 
used to produce stability information useful to ship 
designers and towing operators. Stability information can 
be assembled into a convenient graphical form that clearly 
defines the regions of stable and unstable operation based 
on the parameters the operator has the most control over - 


the placement and length of the towline. 


615, 


For the ship designer, this technique can be useful in 
determining the implications particular design decisions 
would have on the vessel’s performance under tow. 
Depending on the vessel’s use, adjustments to the design 
can be made to improve towing stability, or the customer 
can be forewarned to avoid certain kinds of operations. 
Since nearly all vessels are towed at some time, towing 
performance should be analyzed for all vessels. 

For the towing operator, this technique can provide 
readily available information about how a particular vessel 
will respond under tow. The operator can then adjust the 
towing parameters (i.e., placement the attachment point 
and/or length of the towline) so the tow will be in its 
most stable condition, or, if unavoidable, know that a 
particular towing situation will be potentially dangerous 
and make preparations to deal with it. 

Since ship data is inputted through a data file, the 
towed performance of any vessel can be analyzed with this 
method, including structures such as offshore oil 
platforms. Existing vessels can be analyzed, as well as 


different loading conditions. 
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Two principal disadvantages are associated with this 
technique: 

1. The programs are dependent on the quality of the 
data provided. Determining hydrodynamic coefficients 
and resistance data requires tow tank experiments and 
analysis, and are not obtained for most vessels; 

2. The programs require large amounts of computer time 
and memory to run, which may not be available or too 
costly for potential users, especially to run 
extensive "what-if" scenarios. This problem may be 
alleviated as more inexpensive, high speed, high 
capacity micro- and personal computers become 
available. 

B. RECOMMENDATIONS 

This study was done for only one set of conditions. 
Further research can be done in determining the effect of 
varying conditions, such as different speeds or maneuvering 
by the towing vessel, on the stability of the tow. External 
forces are modelled by the bias in the data file. A 
systematized method of introducing biases into the data 
would enable the analysis of the effect of environmental 
conditions on the towing system. 

Further work should be conducted to improve the "user- 
friendliness" of the programs. As currently configured, the 
programs must be run instructively, and graphics produced 
offline. This is a time consuming process which does not 


use the full capabilities of either the programs or the 


graphics capabilities of the mainframe. Program 
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improvements should focus on streamlining computations and 
user interaction, and incorporating graphics, with the goal 


of making it available as a ship design tool. 


33 


APPENDIX 


Driver programs used in this thesis are shown here. 
Subroutines can be obtained by contacting: 

Prof. F.A. Papoulias 

Department of Mechanical Engineering 

Naval Postgraduate School 

Monterey, CA 93943 


FILE: TOWBIF1 FORTRAN Al 


PROGRAM TOWBI1 TOWOOOIO 

Cc TOWOO00LO 

[ BIFURCATION ANALYSIS OF TOWING SYSTEMS TOWO0020 
Cc PARAMETER DEPENDS ON IPAR TOWO0040 
C IPAR = 1 xP TOWO00S0 

Cc eas 1F TOWO0060 

Cc ats TOWO0070 
Cc IT NEEDS SUBROUTINES FROM TOWING.FTN TOWOOCSO 

es TOWO0990 

IMPLICIT DOUBLE PRECISION (CA-H-0-2) TOWO0100 

DOUPLE PRECISION MASSP-NVD.NV.NRD> ToceNR>»LB>LEN LW, TOWOO110 

] NO.NOU.NOUU.NVVV>NVRR»NVDD NVU,NVUU>NRRR>NRVV> TOWO0120 

2 NRDD.»NRU»NRUU,ND-NDODD»-NDVV.NDRR»NDU,NDUU-NVRD TOWOO120 

Cc TOWO0I40 

DIMENSION IV1(6)+A0656)+VV03) 5X06) WRI6) W106 )+206+6) »SV2(6) TOWO0150 

Cc TOWO0160 

COMMON/INTGR/ISKEG,NREDP,ITYS» ID, IFDS»ISTAB» IPROP TOWOOT70 

COMMON/SPAR/MASSP LW. XPP,YPP.LB TOKOOT8O 

CEMMON/SURGE/SUT?7 } TOWO0190 

COMMON/XSURG/XU»XUU. KXUUU TowocLod 

CCYMON/SWAY/SW(15) TOWOOLTO 
COMMCN/YAW/YA( 16) TOWOOSS 

COMMON/MTER/VCAR -RHO,ABS-CON] .CONC TOWOO0LI30 

COMMON/RESIST/VEL(40) RESI (40) TOWOOLGO 
COMMON/VELE/UEL (100) TOWOOLSO 

COMMCN/POSTHN/X1 6¥1+21 TOWCOLEO 
COMMON/GECM/AL -RW.G,AET.HW- HW] TOWOOL70 
COMMON/PROP/ALE »P,EY-DIA-ANIU TOW90280 
COMMCN/CTNR/XC(99),¥C(99).20(099),TO(99) TOWOO0C30 
COMMON/INTISIC TOWO0300 
COMMON/ DOC /UC. ALPHA TOWOOZI0 
COMMON/UEPT/RLX-RLY>RLO TOWOOZC 
COMMON/SLOPE /PDRXX-PDRXY,PORYX»PORYY TOWO0S2 
COMMON/SLAN/RXX6 »-RYY6+RAX-RYY TOWO0S40 

Cc TOWOOZSO 
OPEN (UNIT=S3S,FILE=*BARGES’ »STATUS="OLD" ) TOWO0260 

OPEN CUNIT=1,FILE="RES0*.STATUS=’NEW® ) TOWO0Z70 

C TOWO0Z80 
OPEN (UNIT=11 -FILE=’RESIR® »STATUS=" KEW" ) TOWD0290 

CPEN (CUNITZIC .FILE=’RESER salella- ten) TOWO0400 

OPEN CUNIT=12.FILE=*RESZR® »STATUS="NEHW’ ) TOWO0410 

OPEN (UNIT=146-FILE=*RESG2’ ,STATUS=’NEW® ) TOWOO0GLO 

OFEN (CUNIT=15-+FILE="RESSR’*® »-STATUS=’NEW® ) TOWO0470 

OPEN CUNIT=16-FILE="RESGER® »STATUS=’NER’ ) TOW00440 
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FOR BIFURCATION ANALYSIS 


»STATUS=' NEW’ ) 
~-STATUS=’ NEW ) 
*STATUS=’ NEW? ) 
~STATUS=’ NEW? ) 
*OTATUS=’* NEW’ ) 
»STATUS=° NEW’ ) 


OF K-TH EQUILIERIUM 


CALL RG(6-6+A.WRIWI-MATZ.2,1V1.SVO,TERI) 


TATION’ »/> 
IATION' ./;> 
TATION’) 


PARAMETER RANGE’) 
NUMBER OF INCREMENTS‘ ) 
EQUILIERIUM NUMEER® ) 


OF EQUILIBRIA® ) 


OPEN CUNT T=2);FILE#"RESII° 
OPEN GUNIT=22 sFICE="RESSI 
OPEN (UNIT T223.FILE=*RESII°’ 
OPEN (UNI T=c4¢.FILE=*°RESGI' 
OPEN GCUNIT=cS-FILE="RESSI° 
OPEN (UNIT=26-FILE=*RESO6I’ 
CALL INPUTD(I0O) 
VCAR =VCAR#1].689D0 
MATZ =0 
IFLOW=! 
WRITE (*,1001) 
READ (4>*) IPAR 
WRITE (51002) 
READ (4%) Al.AS 
WRITE (*,1003) 
READ (4,™) NUMI 
WRITE (*,1005) 
READ (#4) IKB 
WRITE (*.1006) 
READ (™,") NEOL 
IF CIKB.GT.NECL) GO TO 500 
DO 1 I=1.NUM1 
WRITE (=.2001) I.NUM)] 
AAZA]*(AL—Al 8 CIT] 2/7 CNUMI=-1) 
IF CIPAR.EO.1) XPP=AA 
IF (JPAR.EO.2) YPP=AA 
IF CIPAR.EO.3) LW SAA 
AL =LWLBrO.30¢68D0 
ALE=AL 
CALL STABIL(IVV.VV.ISOL) 
SET V=VV(K) 
IF CIVV.NE.NEOL) GO TO 1 
V=VV( IKE) 
Hea CDABS UVIRGT 2.00)! SOR sl) 
CALL EQUILB(V.X-RES-RX-RY) 
CALL LINEAR(X,RES.ARXK-RY) 
DE SCrER ONES 0) STOR ses2 
CALL DEGSTB(DEOS.WR) 
WRITE (1,10) AA.DEOS 
DO 1] J=1.6 
JR=10¢3 
WRITE (JR.10) AA.WROJ) 
<b) (aa oe 
WRITE (41,10) AA WI(J) 
I] CONTINUE 
1 CONTINUE 
500 STOP 
1Q9 FORMAT (2D20.10) 

1001 FORMAT (° ENTER 1 : XP VAR 
l ; 2: YP VAR 
= : 3 : LW VAR 

1092 FORMAT (° ENTER 

1003 FORMAT (° ENTER 

1005 FORMAT (° ENTER 

1006 FORMAT (° ENTER ESTIMATED NO, 

2001 FORMAT (215) 

END 


Spo) 


TCWO0450 
TOWOOG60 
TOWD0470 
TOWO0480 
TCWO0490 
TCWOOSO0 
TCWO0S10 
TOWOOSLCO 
TOWDOS 50 
TCWOOS40 
TOWOOSSO 
TOWOOS 60 
TOWOOS70 
TCWOOSS8O 
TCWOOESO 
TOWO0600 
TOWO0610 
TOWOO6CO 
TOWODGSO 
TCHO06460 
TOWO0ESO 
TOHO066D 
TOWO0670 
TCWO0680 
TCHO0690 
TCWOO700 
TOWo0710 
TCWOO7CO 


TOWOO750 
TOWOO740 
TOWDO750 
TCWO07060 
TCOWOO770 
TCWOO78O 
TOWOO790 
TOWOOSCO 
TCWO0810 
TCWOO8IO 
TCWOOSS0 
TOWOCESO 
TOWO 0850 
TOWOI8ED 
TOWOOS8 70 
TOWO0880 
TOHWO0890 
TCWO0S09 
TOWOO910 
TOWOdsLO 
TOWOOSSO 
TOWOO 940 
TOAOO?S0 
TOWOO 960 
TOWO0970 
TCWOOS80 
TOWNOSIO 
TCwWO1000 
TOWO1010 
TOWO10L0 
TOWO 1020 
TOWO 1040 
TOWO1050 
TOWO1060 
TOWO 1070 
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TOWBIF2 FORTRAN Al 


PROGRAM TOWBI 2 
PROGRAM TOWBIF.FTN 


BIFURCATICN ANALYSIS OF TOWING SYSTEMS 
PARAMETERS ARE: Xp. Yp 
IT NEEDS SUBROUTINES FROM TOWING.FIN 


USER DEPENDENT SUBROUTINES: 
DEGSTB = CURVES ENCLOSING REGION 1I OF FIGLFE I3 
(SUBROUTINE DEGSTB 1S 1N TOWING.FTN) 
Dsl = CURVES ENCLOSING REGION V OF FIGURE 13 
(SUBROUTINE DS1 IS IN SPMBIF) 


IMPLICIT DOUBLE PRECISICN (A-H.O0-2) 

DOUBLE PRECISION MASSP>»NVD»NV>NRD»IZZ>NR»LB-LEN>LW, 
NO.NCU-NOUU.NVVV.NVRR-NVDD-NVU»NVUU,LNRRR>NRVV, 
NRDD»NRU>NRUU,ND»NDDD>NDVV.NDRR  NDU>NDUU,NVRD 


DIMENSION IV1(6),A(6,6)-VV(3)+X(6)-WR(E) WI(6),206-.6),5V2(6) 


COMMON/INTGR/ISKEG,NREDP >» ITYS, ID» IFDS»ISTAB, IPROP 
COMMON/SPAR/MASSP>LW»XPP,YPP LB 
CCMMON/SURGE/SUC7) 
COMMCN/XSURG/KU.XUU, XUUU 
COMMCN/SWAY/SKC1S ) 

COMMON/YAW/YA( 16) 
COMMON/MTER/VCAR + RHO. ABS.CON] »CONZ 
COMMON/RESIST/VEL (40) RESI (G0) 
COMMON/VELE/UEL (100) 
CCMMCN/POSTN/XI5Y1,2} 
COMMON/GEOM/AL.»RW.G.AET» HW. HK] 
COMMON/PROP/ALE.P-.EY-DIA-ANTIU 
COMMCN/CTNR/XC(99),YC(99).20(99),TC(99) 
COMMON/INTI/SIC 

COMMON /DOC/UC -. ALPHA 
COMMON/UEPT/RLX. RLY »RLO 
COMMCN/SLOPE/PDRXX.PDRXY.PDRYX.PDRYY 
COMMON/SLAN/RXX6 -RYYVE RXXRVY 


OPEN (UNIT=SS.FILE=’BSKEGO’ »STATUS=’OLD' ) 
OFEN CUNT T=I1.FILE="’RESIR’ »STATUS='NEW? ) 
OPEN CUNIT=1O.FILE=*RESOR’ ,STATUS=*NEW' ) 
OPEN CUNIT=HIZ,.FILE=*RESIR’ »STATUS="NEW® ) 
OPEN CUNIT=14,.FILE=*RESCR* »STATUS="NEW® ) 
CPEN (UNIT=1S -FILE='’RESSR* »STATUS='NER’ ) 
OPEN CUNIT=16.FILE=’RESoR’* »STATUS='NEW® ) 


CALL INPUTD(10) 

VCAR =VCAR*1.6895D0 
AL sLW*LB¥0.35048D0 
ALE =AL 

MATZ =0 

IFLOW=1] 

EPS @1.0-5 
ILMAX=1500 


WRITE (*,1001) 

READ (*,™) Al,Ao 
WRITE (*,1002) 

READ (*,™) NUM1 

WRITE (*,1002) 

READ (.™) Bl .EB2 
WRITE (*,1004) 

READ (¥.™) NUMO 
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TOWOO010 
TOWOOOZC 
TOWOOOZ6 
ToWo0040 
TOWO00SO 
TOWO0060 
TOWO0070 
TOWO0OZO 
TOWO0090 
TOWO0100 
ToWoo110 
TOWOO1Z0 
TOWO0120 
TOHWO0140 
TOWOO1SO 
TOWOO160 
TOW00170 
TOWO0180 
TOWO0190 
TOWOOZ00 
TOWOOZ10 
Towoor20 
TOWOO2S 

TOWOO24u 
TOHOOZSO 
TOWO0260 
ToWoo270 
TOWOOZEO 
TOWO0S90 
TOWO0Z00 
TOWOOZ10 
TOwooz2 

TOWOOS50 
TOW00240 
TOWO0Z50 
TOHOOZ60 
TOWO0S70 
ToW00380 
TOWODTI0 
TOK00400 
TOWO0410 
TCHOOGZ0 
TOWO0450 
TOW00G40 
TOWO04S0 
TOHOGG60 
TOW00470 
TOWO0480 
TOW00490 
TOWOOSO0O 
TOWOOSIO 
TOWOOS20 
TOWO0S50 
TOKO0S4O 
TOWOOSSO 
TOWO0S60, 
TCWOOS70 
TOW00S80 
TOWO0S 90 
TOWD0600 
TOWO0610 
TOWOCEZ0 
TCHOO63O 
TOWO0e<0 


WRIT 


READ 


WRIT 
READ 
DO 1 
WR 
YP 


XP 


—€ (#,1005) 
(u,a) 1KB 
E (",1006) 
(#,™) IDS 
T=I,NUM1] 
ITE €#,20012 1+NUM] 
Prai*CAc=Alo" (l=) 3/CNUMI-1) 


P=El 


CALL STABILCI1VV,VV,1S0L ) 


SET V=VV(K) FOR BIFURCATION ANALYSIS OF K-TH EQUILIBR1UM 


Vis 


VVC IK8) 


E> CDABSTV).GT.1.D0) STOP. 1111 

CALL EQUILB(V.X-RES-RX RY) 

CALL LINEAR(X,RES-A-RXRY) 

CALL RG(6,6+A-WR WI ,MATO,2,IV1-SV2-IER1 ) 


Wg 
LE 
Le 
DE 
XP 
L 

DO 


(TERT NE .0) STOP cece 
(105.€0.1) CALL DEGSTS(DEOS>WR) 
CIDS.60.<2) CALL DSPCDEOS, KR? 
OSCO=DEO0S 
OO =xXPP 
=0 

2 J=tcosNUMo 
WRITE (",") J 
XPPsBl*(BlI-BljJe(J-1)/(NUMO-1 } 
CALL STABILCIVV>VV>ISOL) 


SET V=VV(K) FOR BIFURCATION ANALYSIS OF K-TH EOQUILIER1UM 


VeVVCIKB) 

IF LPAES(V).GI.1.D0) STOP 111! 
CALL EQUILB(V.X RES RX RY) 

CeLl LINEAR(X.RES+ARXRY) 

CALL RG(6-6-A-WR.W]-MATZ.Z,1V1,SV2-JER] ) 
le SClLen TNE OJ US (Ora ca. 

IF (1D3.€0.1) CALL DEGSTB(DEOS.WR ) 
LE CIES. £0 2) CALL DS PCDEOS uWRD 
DEOSNN=DE0S 

XP NN=XPP 

PR=DEOSO0* DEOSNN 

ip (PR2G1.0.00) GO 10 3 

L=L¢] 

IF (L.GT.6) STOP 1000 

1L=0 

XP O=xXFOO 

XFN=XPNN 

DEOTO=DE0L00 

DEOSN=DEOSNN 

XPL=XPO 

XFR=sXPN 

DEOSL=CPECSO 

DEOCR=CEOSN 

XPP=(XPL*XPR)/2.D0 

CALL STABILCIVV.VV>1SOL) 
V=eVVCIKS) 

TES CDARS(VI<G1..1-00) STOP 1111 
CALL ECUILB(V.X,RES.RX»RY) 

CALL LINEAR(X.RES-ARX-RY) 

CALL RG(6>6+AWROWIMATO.2,1V1,SVO,TERI) 
TP UCIER DWE. 0) “SIOP 22ce 

CALL DEGOTE(DEOS. Wa) 

DEOSM=DEOS 


a7 


TOWO06S0 
TOWO0660 
TOWO0570 
TOW00a80 
TOWO0S90 
TOWOO700 
TOWOO710 
TOWOO720 


TOWOO730 
TOWOO740 
TOWO0750 
TOW00760 
TOWO0770 
TOKO07B0 
TOWOO790 
TOWOOBOO 
TOWOO810 
TOWOOSIO 
TOWNO8SO 
TOHOOBGO 
TOWDO8SO 
TOWO 0860 
TOWO0870 
TOWOO880 
TOWCOB90 
TOWO0900 
TOWOO910 
TCWOO 920 
TOWOO9Z0 
TOWO0940 
TOKOOSSO 
TOWO0960 
TOWCO970 
TOWO0 980 
TOWO09 90 
TOWO1000 
TOWO1010 
TOWO1020 
TOW91020 
TCWO 1040 
TOWO1050 
TOWO1060 
TOWO1070 
TOWO10B0 
TOKO1090 
TOWO1100 
TOWOIJ110 
TOWO1120 
TOWO1130 
TOWO1140 
TCWO1150 
TOWO1160 
TOWO1170 
TOWO1180 
TOWO1190 
TCWO1L00 
TOWO1010 
TOWO L220 
TOWO 1220 
TOWO1240 
TOWO1S50 
TOWD 1260 
TOWO1270 
TOWO1S80 


i) 


10 
100) 
1002 
1003 
1094 
1005 
1006 
200) 


ft) 


XPM=XPP 

PRL=DEOSL *DEOSM 
PRR=DEOSR*DEOSM 

IF C(PRL.GT.0.D0) GO TO 5 
XPO=XPL 

XPN=XPM 

DEOSO=DEOSL 

DEOSN=DEOSM 

TeeIt*1 

1F CIL.GT.ILMAX) STOP 3100 
DIF=DAEBS(XPL- XPM) 

iF (DIF.GT.EPS) GO 70 6 
XP=XPM 

GO TO 4 

1F (PRR.GT.0O.D0) STOP 3200 
XPO=XPM 


XPN=XPR 

DEOSO=CEOSM 

DEOSN=DEOSR 

1E=IEe1 

JF CIL.GT.ILMAX) STOP 3100 
DIF=DABS(XPM-XPR ) 

le (DIF -GT-EPS? GO TO°6 
XP =XPM 

CELaTeee 

WRITE (LLL>10) XP.YPP 
XPOO=KPNN 

DEOSO00=CDEOSNN 


CONT IMLE 
1 CONT 


STOP 
FCR! 
FORM 
FORM 
FORM 
FORM 
FORM 
FORM 
FORM 
END 


SUBR 


IMPLIC1T DOUBLE PRECISION (A-H,0-72) 


DIME 
DEOS 
DO 1 
ie 
DE 
IJ 
CONT 
DEOS 
DO 2 
IF 
Ie 
DE 
IJ 
CONT 
DEOS 
DO = 
1F 
Miss 
CONT 


INUE 


AT (2020.10) 


AT (° ENTER RANGE OF Yp VARIATION’ ) 


AT (° ENTER NUMBER OF INCREMENTS 


INS Vip.) 


AT (* ENTER RANGE OF Xp VARIATION® ) 
AT (* ENTER NUMBER OF INCREMENTS IN Xp° ) 


AT (° ENTER EQUILIBRIUM NUMEER' ) 
AT (* ENTER DEGREE OF STABILITY CONTROL" ) 
‘AT (215) 


OUTINE DSI] (DEOS.WR) 


NSION WR6) 
1=-1.D30 
121.6 
(WRO1).LT.DEOS]) GO TO } 
OS1]=WROI) 
=] 
INVE 
22-1.D30 
T=1,6 
C(1J.£0.1) GO TO 2 
CHREI).LT.DEOSS) GO TO 2 
OS2=hwR(1) 
J=] 
INUE 
z-].050 
T=1+6 
CieEG. 1 J0R.1.EO.1 JJ) a0 


(WR(13.GE.DEOS) DEOS=WR(1) 


INVE 


RETURN 


ED 
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TOWO1290 
TOWO1200 
TOWO1210 
TOWO132 

TOWO132 

TCWO1 340 
TOWO1250 
TOWO1260 
TOWO1370 
TOWO1I80 
TOWO1390 
TOWO1400 
TOWO1410 
TOWO1420 
TOWO14350 
TOWO1440 


TOWO1450 
TOWO 1460 
TOW01470 
TOWO1480 
TOWO1490 
TOWO1500 
TOWO1510 
TOWO1520 
TOWO1550 
TOWO1540 
TOWO1550 
TOWO1560 
TOWO1570 
TOWOI1ES80 
TOKO1S90 
TOWO1600 
TOWO1610 
TOWO1ECO 
TOWO1550 
TGWO01640 
TOWO1050 
TOWO1660 
TOWO1E70 
TOWO1680 
TOWO1690 
TOWO1700 
TOWO)710 
TOWOL?Z<2 

TOWO1750 
TOWO1 740 
TOWO1750 
TOWO) 760 
TOWO1 770 
TOWOI780 
TOWO1790 
TOWO1800 
TOWOI810 
TOWO1820 
TOWO1820 
TOWO1840 
TOWO1850 
TOWO1860 
TOWO1870 
TOWO1880 
TOWO1890 
TOWO1900 
TCWO1910 
TOWO) 820 


FILE: 


ono aananqananaqaaqagnn 


TOWBIFS 


PROGRAM 


FORTRAN Al 


TOWSIF3.FT™N 


BIFURCATION ANALYSIS OF TOWING SYSTEMS 
PARAMETERS ARE: Xp> Lw 
IT NEEDS SUBROUTINES FROM TOWING.FIN 


USER DEPENDENT SUBROUTINES: 


DEGSTE 


CURVES ENCLOSING REGICN II OF FIGURE 13 
(SUEROUTINE DEGSTB IS IN TOWLNG.FIN) 
CURVES ENCLOSING REGION V OF FIGURE 13 
(SUBROUTINE DS! IS IN SPMBIF) 


IMPLICIT DOUBLE PRECISION (A-H,0-2) 


DOUBLE PRECISION MASSP,NVD.NV,NRD,IZ2,NR>LB-LEN-LW, 
NO,NOU,NOUU.NVVV-NVRR»>NVDD>NVU.NVUU-NRRR>NRVV, 
NRDD,»NRU,NRUU,ND.NOOD»NDVV,NDRR»NDU.-NDUU,NVRD 


DIMENSION 1V1(6)5A(6.6)-VV(3).X(6)»>WR(6)W1(6).2(656)2+SV2(6) 


COMMON/ INTGR/ISKEG.NREDP,ITYS.1D>IFDS,ISTAB, IPROP 
CCHMON/SPAR/MASSP »>LW, XPP>YPP,LB 

COMMON/ SURGE /SU(7) 
COMMON/XSURG/XU> XUU» XUUU 

COMMONS SWAY/SWO15 ) 

COMMON/YAW/YA( 16) 

COMMON/MTER/VCAR »RHO,ABS,CON] ,COND 
CCMMON/RESIST/VEL(40),RESI(G0) 
COMMON/VELE/UEL (100) 
COMMON/POSTN/XI-Y1>21 

CCMMON/GEOM/AL .RW,G,AET, HW, HW) 
CCMMON/PRCP/ALE>P.EY-DIA,ANLU 

COMMON /CINR/AC( 89), YCC99 7, 2C(99)5TCC9S) 
COMMONS INT1/1C 

CCMMON/DOC/UC »ALPHA 
CCMMON/UEPT/RLX-RLY-RLZ 
COMMCN/SLOPE/PDRXX»+PDRXY,»PDRYX>»PDRYY 
COMMON/SLAN/RXX6 »RYYO>RKXXRVY 


OPEN (UNI1T=35S,FILE=’ TANKERO* ,STATUS='OLD’ ) 
CPEN CUNIT=1.F1LE="’RESIR' »STATUS=°* NEW’ ) 
OPEN (UNIT=2,FILE=’RESOR' »STATUS=’NEW’ ) 
OPEN (UNI1T=3.FILE="’RES3R’ »STATUS=’ NEW’ ) 
OPEN (UNI T=4,FILE="RESGR'»STATUS="NEW® ) 


CALL INPUTD(10) 
VCAR =VCAR«1.689D0 
AL  sLW*LB*0.3068D0 


ALE =AL 
MATZ =0 
IFLOW=1 


ILMAX=1500 


EPS 21. D=5 


WRITE (9.1001) 


READ (m,*) 


Al.A2 


WRITE (*",1002) 


READ (".*) 


NUM} 


WRITE (*,1002) 


READ ("5") 


Bl.E2 


WRITE ("+1004 ) 


READ (#545) 


NUMS 


WRITE (",1005) 


READ (#.™} 


IKB 


ao 


TOWO0610 
TowO00ZD 
TCWOO0Z0 
TOWO0040 
TOWO00S0 
TOWOOD6O 
TOWO0070 
TOW00080 
TOWO0090 
TOWO0100 
TOWCO11O 
TOWOO1CO 
TOWO0130 
TOWO0140 
TOWoo15O 
TOWO0160 
TCWOO170 
TOWOO!80 
TOWOO1290 
TOWOOZ00 
TOWOO210 
Towoos 

TOWOD23 

TOWODS40 
TCWOO2S0 
TCWO0=60 
TOWO0270 
TCHOOZO 
TCWOOZ90 
TOW00200 
TCWO0310 
TOWOOSZ0 
TOWO033 

TOWO0340 
TOWOOD350 
TOWO0T60 
TCWOOZ70 
TOWCO380 
TCWO0S90 
TOWO0600 
TOWOOG10 
TOWOO420 
TOWO0S50 
TOWO0440 
TOwo04S0 
TOWD0460 
TOWONS70 
TOWOO<80 
TOWO0490 
TOWOOS00 
TOWOOS10 
TOWOOS=0 
TOWOOS50 
TOWOOS<O 
TOWOCESO 
TOWOOS6O 
TOWOCS70 
TOWOOSS0 
TOWOCS°O 
TOWC0600 
TOWO0610 
TOWOCELO 
TOWOOSSO 
TOWD0660 


oO 


” 


40 


WRITE (*.1006) TOWOCESO 
READ (m.™) 1DS TOWO0660 
DO 1 1l=1].NUM!] TOWDD670 
WRITE (4.2001) I>NUM1 TOWOD6RBO 
LW =Al*(AD-Al)*(I-1)/(0NUMI-1) TOWD0690 
AL =LW*L80.3048D0 TOW00700 
ALE=AL TOWO0710 
XPP=B1 Towoo720 
CALL STABILCIVV,VV,1SOL) TOWD0730 
TOWO07460 
SET V=VV(K) FOR BIFURCATION ANALYSIS OF K-TH EOQUILI8R1UM TOWO0750 
TOW00760 
V=eVVCIKB) TOWO0770 
IF (DABS(V).GT.1.D0) STGP 1122 TOWO0780 
CALL EQUILB(V-.X-.RES-RX.RY) TOW00790 
CALL LINEAR(X,RES+AsRX-RY) TOWOO8O0 
CALL RG(6-6-A-WRIWI-MATI+2+1V1»SV2>1ER1 } TOWO0810 
IF CIER1.NE=0) STOP e2cec TOWOOB=O 
IF (I1DS.EO.1) CALL DEGST8(DECS.WR) TOWOO08Z0 
IF (1DS.E0.2)} CALL DS1(DEOS.bR) TOWO0840 
DEOSO0=DE0S TOWOOSBSO 
XPOO =XPP TOWO0860 
iL =0 TOWO0870 
DO 2 J=2,NUM2 TOWOOSSO 
XPP=B1*( 82-81 )"(J-1)/(NUMO-1 ) TOWOC890 
CALL STABILCIVV,VV,ISOL) TOWO0900 
TOWO0910 
SET V=VV(K) FOR SIFURCATICN ANALYSIS OF K-TH EQUILIBRIUM TOWOON2920 
TOWO09S0 
VeVVC IKE) TOWND940 
LE SCORESUV I 2G in DOs oer s TOWO0950 
CALL EQUILB(V.X-RESRXRY) TCW90900 
CALL LINEAR(X.RES»A.RX-RY) TOWO0970 
CALL RG(6-6-A-WR WI .MATZ.2-1V1,SV2-IERI ) TOWD0 S80 
1F (TER) .NE.0) STOP Gece TOWN0990 
TES GILDS SEC CAE DEGSTIEIGEeOS..NRD TOWO1000 
LF CIOS EO lS) CALL DSC pegs HAD TOWO1010 
DEOSNN=DEOS TCWO1020 
XPNN=XPP TOWO1039 
PR=DEOTOO*DEOSNN TOWNO10<60 
LE SCRREGT 0. DOD: GOO: TOWO1050 
Cel] TOWO1060 
IF (L.GT.4) STOP 1000 TCWO1070 
IL=0 TOWO10380 
XPQ=XPOUO TOWO1090 
XPN=XPNN TOWO1100 
CEOLO=CEOSO0O TOWO1110 
DECSN=DEOSNN TOWO1120 
XPL=XPO TOWO1130 
XFR=XPN TOWO1140 
DEOSL=DEOSO TOWO1150 
DEOSR=DEOSN TOWO1160 
XPP=(XPL*XPRI/2.D0 TOWO1170 
CALL STABILCIVV-VV.ISOL) TOWO1180 
VeVVCIKE) TOWO11°9O 
LF ADASSUV).G1.1.009) SeOr Tit) TOWOI200 
CALL EQUILS(V-X.RES-.RX.-RY) TOWO1210 
CALL LINEAR(X-RES-A+RX,RY) TOWO1220 
CALL RG(6-6,A,WR.W] MATS, 2+1V1,SV2,. LER] ) TOWO123 
IF ClER] NE. 03° Si1Gr cece TOWG1240 
CALL DEGST8(DEOS-kR) TOWO1250 
DEOSM=DECS TOWOIZED 
XPM=XPP TOWO1270 


1002 
109¢ 
1005 
1004 
1005 
1006 
200} 


PRL=DEOSL*DEOSM 
PRR=DEOSR*CEOSM 

IF (PRL.GT.0.D0) GO TO 5 
XPO=XPL 
XPN=XPM 
DEOSO=DEOSL 


DEOSN=DEOSM 


IL=IL*? 
IF CIL.GT.ILMAX) STCP 3100 
DIF =DABS (XPL-xKPM) 
JESCOie GT EPS) GOTO 6 


XP =XPM 


GO TO 4 
IF CPRR.GT.0.D0) STCP 3200 
XPO=XPM 
XPN=XPR 
DEOSO=DE0SM 


DEOSN=DEOSR 
Me=lEet 

IF CIL.GT.ILMAX) STOP 3200 
DIF=DABS(KPM—XPR ) 

IF (DIF.GT.EPS) GO TO 6 


XP =XPM 


WRITE (L~10) XP>LW 
xPOO=KPNN 
DEOSO0=DEOSNN 


CONTINUE 


CONTINUE 
STOP 


FORMAT (2020.10) 


FORMAT (° 
FOPMAT (° 
FORMAT (° 
FORMAT (° 
FORMAT (’ 
FORMAT {° 


ENTER 
ENTER 
ENTER 
ENTER 
ENTER 
ENTER 


FORMAT (2135) 


END 


RANGE OF Lw VARIATICN®* ) 
NUMEER OF INCREMENTS IN Lw') 
RANGE OF Xp VARIATION’ ) 
NUMEER OF INCREMENTS IN Xp‘) 
EQUILIERIUM NUMBER’ ) 

DEGRES OF STABILITY CONTROL") 


SUBROUTINE DSI (DEOS-bR) 

IMPLICIT DOUBLE PRECISION (A-H-O-¢) 
DIMENSION WR(6) 
DEOS)}2-1.Dz0 


DO 1 121-6 


IF (WRC(I).LT.DEOS!) GO TO 1 
DEOS1=WRC1) 


IJ=! 
CONTINUE 


DEOS2 2-1 .030 


pO 2 J=1>6 


TE CIS.EO. 1) GO_10 < 
1] CUP C1). 1. Decs. 3 GO-TO: < 
DEOSZ=WR(1) 


IJJj=1 
CONTINUE 


DEOS=-1.D=0 


DO 2 1=1>+6 


1F €1.€0.13.0R.1.£9.iJJ) GO TO 5 
IF (WR(1).GE.DECS) DEOS=WR(1) 


CONTINUE 
RETURN 
EbS 


4) 


TOWOI 280 
TOWO 1290 
TOWO 1500 
TCWOII10 
TCWO S20 
TOWOI3: 

TOWO1540 
TOWO1S50 
TOWO1260 
TOWO1570 
TOWO1380 
TOWO15390 
TOWO1400 
TOWO1410 
TOWO1420 
TOWO14G30 
TOWO1440 


TOWO1450 
TOWO1460 
TOWO1470 
TCWO1489 
TOWO14°0 
TOWO}500 
TOWO1510 
TOWO1S20 
TOWOIS20 
TOWO1540 
TOWOIESO 
TOWOI550 
TOWO1570 
TCWO1580 
TOWO 1590 
TOWO16C0 
TCWO161C 
TOWOL 600 
TOWO1630 
TOKO} 640 
TOWO1650 
TOWO1660 
TOWO1679 
TOWO16E0 
TCWO1690 
TOWOL7 CO 
TOWO1710 
TOWO17=20 
TOWOL730 
TCWO1L760 
Towo1750 
TCWO1760 
TOWO1770 
TOWOL FEO 
TOWD1 790 
TOWD1800 
TCWO1810 
TOWO1820 
TOWO1S50 
TOWD1840 
TOWO1850 
TOWD18 069 
TOWO1870 
TCWOI8S8d 
TOKO1890 
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